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I.  Introduction 

.  .v 

'  A  lot  of  work  has  been  done  towards  the  direction  of  developing  communication 
protocols  that  determine  how  a  single  common  resource  can  be  efficiently  shared  by  a 
large  population  of  users.  By  now,  it  is  well  known  that  fixed  assignment  techniques  are 
not  appropriate  for  a  system  with  large  population  of  bursty  users.  In  the  latter  case,  ran¬ 
dom  access  protocols  are  more  efficient  and  many  of  them  have  been  suggested  [1],  [2], 
Usually,  the  amount  of  information  transmitted  per  time  is  of  fixed  length,  called  a 
packet.  In  most  of  the  systems,  time  is  divided  into  slots  of  length  equal  to  the  time 
needed  for  a  packet  transmission  (slotted  systems). 

The  deployment  of  an  ever  increasing  number  of  multi-user  random  access  com¬ 
munication  networks,  brought  up  the  question  of  how  packets  whose  destination  is 
another  network,  should  be  handled.  Thus,  the  issue  of  network  interconnection  or 
multi-hop  packet  transmission,  ariseS|[3],  [6],  [7], 

The  basic  problem  in  analyzing  interconnected  systems  is  that  of  characterizing  the 
output  process  of  a  multi-user  random  access  communication  system;  i.e.,  the  departure 
process  of  the  successfully  transmitted  packets.  Another  problem  is  how  a  random 
access  protocol  operates  in  the  presence  of  a  node  that  forwards  exogenous  traffic  coming 
from  other  networks.  The  latter  problem  can  be  avoided  by  assigning  a  separate  channel 
to  the  exogenous  traffic.  In  this  case,  the  operation  of  the  system  is  not  affected  by  the 
exogenous  traffic  but  the  problem  of  optimum  allocation  of  the  available  resources  (chan¬ 
nels),  arises.  The  latter  issue  has  been  discussed  in  [3],  where  the  objective  is  to  maxim¬ 
ize  the  throughput  of  the  interconnected  networks.  In  [3],  delay  analysis  was  not  per¬ 
formed  and  only  simulation  results  were  obtained. 


The  output  process  of  a  multi-user  random  access  communication  system  depends 
on  the  protocol  that  has  been  deployed.  Description  of  that  process  is  a  difficult  task  and 
only  approximations  based  on  special  assumptions  have  been  attempted,  [4]-f7], 

In  this  paper,  the  output  process  of  a  multi-user  random  access  slotted  communica¬ 
tion  system  is  approximated  by  a  process  which  depends  on  an  underlying  3-state  Mar¬ 
kov  chain  and  a  stationary  mapping  rule.  A  visit  to  a  state  takes  place  at  the  end  of  a  slot 
and  corresponds  to  one  or  zero  departures,  depending  on  the  state  and  according  to  the 
stationary  mapping.  The  state  space  of  the  Markov  chain  is  given  by 

S  =  {0,  1,  2}  =  {I,  S,  C} 

where  1(0),  S(l)  and  C(2),  respectively,  correspond  to  an  idle,  involved  in  a  successful 
packet  transmission  or  involved  in  a  packet  collision,  slot.  The  feedback  information 
available  to  the  users  at  the  end  of  each  slot  may  not  be  ternary. 

The  motivation  behind  the  adoption  of  the  Markov  model  to  approximate  the  output 
process  is  due  to  the  fact  that  the  output  of  the  channel  at  the  current  time  slot  depends 
strongly  on  the  outcome  in  the  previous  slot.  This  dependence  is  introduced  by  the  colli¬ 
sion  resolution  algorithm  that  is  always  present  in  a  random  access  multi-user  communi¬ 
cation  system.  Thus,  an  independent  packet  departure  model  would  be  both  unrealistic 
and  intuitively  not  pleasing.  On  the  other  hand,  a  2-state  (arrival-no  arrival)  Markov 
model  would  also  be  inefficient  since  idle  and  collided  slots  would  correspond  to  the 
same  state.  Clearly,  the  behavior  of  the  system  in  a  slot  following  a  no-arrival  slot 
depends  strongly  on  whether  the  no-arrival  slot  was  an  idle  or  a  collided  one.  The  latter 
observation  motivates  the  selection  of  the  proposed  3-state  Markov  approximation  of  the 


output  process. 


In  the  next  Section,  the  approximation  model  for  the  output  process  is  introduced. 
The  algorithm,  which  is  used  as  an  example  in  this  paper  for  conflict  resolution  within 
the  network,  is  also  described  in  that  section. 

In  section  HI,  the  steady  state  and  state  transition  probabilities  of  the  underlying 
Markov  chain  are  calculated  as  the  average  state  and  state  transition  probabilities  of  the 
channel  status,  under  the  operation  of  the  deployed  random  access  algorithm. 

In  section  IV,  the  queueing  problem  that  appears  in  a  central  node  which  receives 
and  retransmits  packets  originating  from  N  independent  multi-user  random  access  slotted 
communication  systems,  is  formulated.  Then,  a  method  is  provided  to  calculate  the 
mean  time  that  a  packet  spends  in  the  central  node.  The  latter  is  used  as  a  measure  of  the 
performance  of  the  proposed  approximation  model  of  the  output  process. 

In  the  last  section,  results  for  the  average  time  that  a  packet  spends  in  the  central 
node  are  presented  and  conclusions  are  drawn.  The  performance  of  the  proposed  approx¬ 
imation  of  the  output  process  is  compared  with  that  of  a  2-state  Markov  approximation 
and  that  of  an  independent  packet  departure  process.  Simulation  results  of  the  actual  sys¬ 
tem  are  also  provided. 

II.  The  approximation  of  the  output  process 


Let  {x\  }j>0  denote  the  process  that  describes  the  state  of  the  ith  channel  at  the  end  of 
the  slots.  The  state  space  of  [x]  }j>o  is  assumed  to  be  given  by 


S‘  =  {0,  1,2}  =  {I,  S,C} 


(1) 


as  it  is  explained  in  the  Introduction.  Since  a  packet  appears  in  the  output  process  if  and 
only  if  it  is  the  only  one  transmitted  within  the  corresponding  slot,  the  output  process  of 
the  i1*1  network,  [a]  }j>0,  can  be  clearly  described  via  the  mapping 


1  if  x]  =  1 
0  if  x]  =  0,  2 . 


(2) 


The  process  {xj }  j>0  is  controlled  by  the  deployed  random  access  algorithm  and  it  is 
generally  non-Markov.  However,  this  process  can  be  approximated  by  a  Markov  process 
{x]}j>0  which  has  the  same  state  space  as  {x]}j>o  and  is  ergodic  within  the  stability 
region  of  the  random  access  algorithm.  A  Bernoulli  model  of  the  output  process  is 
unrealistic  due  to  the  dependence  introduced  by  the  collision  resolution  algorithm.  A  2- 
state  Markov  model  can  also  be  inefficient  for  the  reasons  mentioned  in  the  Introduction. 


As  an  example,  we  consider  multi-user  random  access  slotted  communication  net¬ 
works  in  which  a  binary  (Collision/Non-Collision)  feedback  limited  sensing  collision 
resolution  algorithm  is  deployed.  This  algorithm  has  been  developed  and  analyzed  in 
[11]  and  [10].  There,  analysis  was  limited  to  the  derivation  of  the  maximum  stable 
throughput  and  the  average  packet  delay.  The  characterization  of  the  process  of  the  suc¬ 
cessfully  transmitted  packets,  i.e.  the  output  process  of  the  network,  is  still  an  open  prob¬ 
lem. 


A  brief  description  of  the  collision  resolution  algorithm  is  provided  at  this  point. 
Each  user  is  assigned  a  counter  whose  initial  value  is  zero  (no  packet  to  be  transmitted). 


This  counter  is  updated  according  to  the  steps  of  the  algorithm  and  the  feedback  from  the 
channel.  Upon  packet  arrival,  the  counter  content  increases  to  one.  Users  whose  counter 
content  is  equal  to  one  at  the  beginning  of  a  slot,  transmit  in  that  slot.  If  the  channel 
feedback  is  collision  (C),  the  counters  whose  content  is  greater  than  one  increase  it  by 
one;  the  counters  whose  content  is  one  maintain  this  value  with  probability  p  (splitting 
probability)  or  increase  it  to  two  with  probability  1-p.  If  the  channel  feedback  is  non¬ 
collision  (NC),  all  non-zero  counters  decrease  their  content  by  one.  A  detailed  descrip¬ 
tion  of  the  algorithm  can  be  found  in  [10],  [1 1]. 

To  be  able  to  apply  the  proposed  3-state  Markov  approximation  of  the  output  pro¬ 
cess  of  a  network  that  operates  under  the  described  collision  resolution  algorithm,  the 
steady  state  and  state  transition  probabilities  of  the  output  process  need  to  be  calculated. 


III.  The  steady  state  and  state  transition  probabilities 


Since  the  Markov  model  is  only  an  approximation  of  the  output  process,  it  seems 
natural  to  estimate  the  steady  state  and  state  transition  probabilities  of  the  Markov  chain 
by  calculating  the  steady  state  probabilities  that  a  particular  state  or  state  transition 
occurs  in  the  output  process,  under  stable  operation  of  the  network.  This  calculation  is 
not  always  straightforward.  The  procedure  to  be  followed  depends  on  the  class  of  ran¬ 
dom  access  algorithms  which  has  been  employed. 

In  this  section,  we  calculate  these  probabilities  for  a  network  in  which  a  binary 
(C/NC)  feedback  limited  sensing  random  access  algorithm  has  been  deployed.  The  pro¬ 
cedure  that  is  followed  can  also  be  applied  to  a  system  in  which  some  other  limited  sens¬ 
ing  algorithm  is  deployed.  In  the  latter  case,  the  complexity  of  the  calculations  may  be 


increased. 


An  important  quantity  in  the  analysis  of  the  random  access  algorithm  under  con¬ 
sideration,  is  the  session  length.  A  technical  definition  of  the  session  via  the  use  of  an 
imaginary  marker,  can  be  found  in  [12],  [13].  In  essence,  a  session  is  a  sequence  of  con¬ 
secutive  slots  that  starts  and  ends  at  points  in  which  the  system  regenerates  itself.  The 
multiplicity  of  a  session  is  determined  by  the  number  of  packet  transmission  attempts  at 
the  first  slot  of  a  session. 

The  following  quantities  are  useful  in  the  analysis  that  is  presented  in  this  section, 
ij-slot :  A  slot  that  is  in  state  i  and  it  is  followed  by  a  slot  in  state  j,  i,  j  e{I,  S,  C}. 

i-slot :  A  slot  that  is  in  state  i,  i  £{I,  S,  C}. 

internal  slot  :  An  ij-slot  is  internal  if  both  slots  belong  to  the  same  session;  i-slots  are 
always  internal,  i,  j  £  {I,  S,  C}. 

lk  :  Length  of  a  session  of  multiplicity  k  (in  slots). 

Lk  :  Expected  value  of  lk. 

L  :  Expected  value  of  Lk  with  respect  to  k. 

x]J  (XL)  :  Number  of  internal  ij-slots  (i-slots)  in  a  session  of  multiplicity  k, 

i,  j  £  {I,  S,  C}. 

T[!  (Tk) :  Expected  value  of  xj|  (xk). 

T'j  (T1)  :  Expected  value  of  Tj^  (Tk)  with  respect  to  k. 

t„J  On)  :  Number  of  internal  ij-slots  (i-slots)  in  the  nth  session  from  the  time  origin. 

tn  :  Length  of  the  nIh  session  from  the  time  origin  (in  slots). 
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q1  (q'j) :  Steady  state  probability  of  a  slot  in  state  i  (steady  state  probability  of  a  slot 

in  state  i  followed  by  a  slot  in  state  j),  i,  j  e  (I,  S,  C}  . 

p(i,j):  Conditional  steady  state  transition  probability  of  having  a  j-slot  in  the  next 

slot  given  that  we  have  an  i-slot  in  the  current  slot,  i,  j  E  {I,  S,  C}. 

The  objective  in  this  section  is  the  calculation  of  the  conditional  and  state  probabili¬ 
ties  p(i,j)  and  q1  ,  respectively,  i,  j  e  {I,  S,  C}.  Clearly,  p(i,j)  can  be  obtained  by  divid¬ 
ing  the  joint  probability  q1-*  by  the  state  probability  q1,  i,  j  £  (I,  S,  C).  Thus,  it  suffices 
to  calculate  the  probabilties  q1*  q1-* ,  i,  j  £  {I,  S,  C}. 

Under  stable  operation  of  the  algorithm,  the  steady  state  probability  of  having  a  slot 
involved  in  a  successful  transmission,  qs,  in  the  output  process,  is  simply  given  by  the 
cumulative  input  rate  to  the  network,  X.j.  It  should  be  noted  that  the  validity  of  the  results 
that  are  obtained  in  this  section  is  contingent  upon  the  validity  of  the  assumption  that 
A.j  <  Sif  max-  where  Sj_  max  is  the  maximum  stable  throughput  of  the  algorithm  associated 
with  the  i*  network. 

As  it  will  become  clear  later,  an  important  quantity  for  the  calculation  of  the  proba¬ 
bilities  q'.q'J,  i,j  e  (I,  S,  C},  is  the  mean  session  length,  L.  The  latter  can  be  calculated 
by  following  procedures  similar  to  those  that  appear  in  [1 1],  [12],  [13],  [14].  In  fact,  for 
the  specific  algorithm  under  consideration,  L  has  been  calculated  in  [10]  and  [11],  and 
numerical  results  can  be  found  in  [10].  For  the  completion  of  the  calculations  of  this  sec¬ 
tion  we  start  by  calculating  L. 

From  the  description  of  the  algorithm  the  following  equations  can  be  written. 


lo  =  *  ,  Ii=l 


(3a) 


Fi  and  F2  are  Poisson  random  variables  over  T=1  (length  of  a  slot)  with  probability  den¬ 
sity  Pf( )  ;  4>i  is  Binomial  with  parameters  k  and  p  (p  =  .5)  and  probability  density  bk( ) . 
By  considering  the  expected  values  in  (3)  with  respect  to  all  random  variables  involved, 
we  obtain 

L0  =  1  ,  Lj  =  1  (4a) 

Lk  =  1  +  X  X  Pf(Fi)  bk  (4>i)  L<d,+f,  +  X  X  Pf(F2)  bk  (0i)Lk-o1+F2  *  k>2.  (4b) 

F,=0  6,=0  F2=0<>1=0 

The  infinite  dimensionality  linear  system  of  equations  in  (4)  can  be  written  in  the  general 
form 

Lk=hk  +  £akjLj  ,  k>0.  (5) 

j=o 

The  most  widely  used  definition  of  stability  is  the  one  which  relates  it  with  the 
finiteness  of  Lk,  for  k<°<>.  In  [10],  [11]  it  has  been  found  that  the  system  is  stable  for 
ki<Si,  max  -  -36  (packets/packet  length).  The  authors  in  [10],  [11]  were  actually  able  to 
find  a  (linear)  upper  bound  on  Lk  which  is  finite  for  k<°°.  max  is  then  defined  as  the 
supremum  over  all  rates  X*  for  which  such  a  bound,  Lk,  was  possible  to  obtain. 

The  existence  of  Lk  <°°,  for  k<°o,  implies  that  (5)  has  a  non-negative  solution,  Lk; 
the  solution  of  the  finite  dimentionality  system  of  equations 

J 

Lk  =  hk  +  X  akj  lj  .  (6) 

j=o 

is  a  lower  bound  on  Lk  and  Lk  — >Lk  as  J— >«»,  ( 1 1 1,  [  14],  [15], 

It  turns  out  that  for  sufficiently  large  J  (e.g.  15),  Lk  is  extremely  close  to  Lk  and  thus, 
for  practical  purposes,  Lk  is  considered  to  be  equal  to  Lk,  especially  for  kt  outside  the 
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neighborhood  of  Sj<  max  .  (The  latter  can  be  shown  by  calculating  a  tight  upper  bound  on 
Lk  and  observing  that  it  almost  coincides  with  Lk,  see  [1 1],  [13],  [14]  for  the  procedure). 

By  solving  (6)  with  akj  and  hk  given  by, 

aoj  =  atj  =  0  ,  0<j<J  (7a) 

min{j,k}  k 

akj  =  Z  pf0~4>i)  t>k  (<t>i)+  X  Pf0“k-H>i)  bk  (4>i)  ,  0<j<J  (7b) 

0,=O  <t>!  =  max{k-j,0} 

and 

hk  =  1  ,  0<k<J  (8) 

we  calculate  the  mean  session  length  of  multiplicity  k.  Since  the  multiplicities  of  succes¬ 
sive  sessions  are  independent  and  identically  distributed  random  variables,  the  mean  ses¬ 
sion  length,  L,  is  calculated  by  averaging  Lk  over  all  k  ;  k  is  the  number  of  arrivals  in  a 
slot  from  a  Poisson  process  with  intensity  Xj.  In  fact,  the  average  for  k<J  is  sufficient. 

In  the  next  three  sub-sections  we  calculate  the  steady  state,  q\  and  joint,  q1J,  proba¬ 
bilities,  i,  j  e  (I,  S,  C). 

Ill- A:  Calculation  of  q1,  qCI,  qcs. 

We  start  by  calculating  the  joint  probability  qCI;  i.e,  the  probability  that  a  slot  is 
involved  in  a  packet  collision  and  the  slot  that  follows  is  idle.  By  using  the  notation 
which  was  introduced  at  the  beginning  of  this  section,  the  following  equations  can  be 
written. 

Xo‘=0  ,  if  =0  (9a) 

xk'  =  1{0,+F,=0}  +  X0,I+F1  +  Tk-H5,+F2  >  k>2.  (9b) 

The  definition  of  the  session  implies  that  the  last  slot  of  a  session  cannot  be  involved  in  a 


collision.  Thus,  all  Cl-slots  are  internal  slots  and  the  equations  in  (9)  take  into  considera¬ 
tion  all  possible  Cl-slots  that  may  appear  in  the  output  process. 

By  considering  the  expected  values  in  (9)  with  respect  to  all  random  variables 
involved,  we  obtain  an  infinite  dimensionality  linear  system  of  equations  with  respect  to 
Tk'  which  is  of  the  form  of  the  one  in  (5).  A  truncated,  up  to  J,  version  of  that  system 
would  be  of  the  form  of  that  in  (6).  In  fact,  the  coefficients  akj ,  0<k<J ,  0<j<J,  are  also 
given  by  (7),  while  the  constants  hk  ,  0<k<J  are  given  by 

hg1  =  hf  =  0  ,  h£‘  =  Pf(0)  bk  (0)  ,  2<k<J.  (10) 

Since 


Tg‘<Lk  (11) 

and  since  Lk<°°,  for  k<°°  and  inside  the  stability  region  of  the  algorithm,  it  is  implied 

that  the  finite  system  of  equations  with  respect  to  Tk(,  0<k<J,  has  a  non-negative  solu- 

tion  Tk  ,  0<k<J,  which  converges  to  Tk  as  J— >°°  The  comments  about  the  rate  of  the 
“  “Cl 

convergence  of  Lk  are  valid  for  Tk  also.  Thus,  an  extremely  good  approximation  of  the 
mean  number  of  Cl-slots  in  a  session,  T"1,  is  obtained  by  averaging  Tk  over  all  0<k<J. 

By  applying  the  law  of  large  numbers  we  obtain  the  following  expression  for  the 
joint  probability  qCI. 


qC,=  lim 


M 

I  tg 

n=l 


Cl 


M— 


M 


Itn 


=  lim 
M 


±  y  ,CI 


1 


M 


(12) 


M  £ 


(^1 

Clearly,  the  random  variables  t„  ,  n>l,  are  independent  and  identically  distributed.  The 


same  holds  for  tn  ,  n>  1 .  Since 


12 


! 


> 


i 


E{t£I}=TCI  <oo  and  E{tn}  =  L<«» 


(13) 


by  applying  the  strong  law  of  large  numbers  to  (12)  we  obtain,  [18], 


ci  T° 

q  =T 


(14) 


To  calculate  the  probabilities  q1  and  qcs  we  follow  exactly  the  same  procedure. 
and  Tjp  are  computed  by  solving  finite  dimensionality  linear  systems  of  equations  of  the 
form  of  (6),  whose  coefficients  are  derived  in  Appendix  A.  Then,  T1  and  T05  are  com¬ 
puted  as  in  the  case  of  T01.  Finally,  the  strong  law  of  large  numbers  is  again  applied  and 
i  r\ 

the  probabilities  q  and  q  can  be  calculated  by  the  following  expressions. 

T1  rs  T055 


qi=J_  qcs  _ 

4  L  ’  4  L 


(15) 


There  is  an  easier  way  to  calculate  the  steady  state  probability  of  having  an  idle  slot, 
q1,  for  the  particular  algorithm.  This  procedure  utilizes  a  property  of  the  binary  rooted 
tree  which  can  describe  the  splitting  operation  of  the  algorithm.  According  to  that  pro¬ 
perty,  the  number  of  terminal  nodes  is  equal  to  the  number  of  intermediate  nodes  plus 
one.  Since  the  intermediate  nodes  correspond  to  collisions  and  the  terminal  nodes 
correspond  to  idle  or  successful  slots,  within  a  session,  the  following  equation  holds. 


E+S- 1 =C 

where  E,  S,  C  is  the  mean  number  of  idle,  successful  or  collided  slots  in  a  session.  Since 
S  =  XjL  and 


E  +  S  +  C  =  L 

we  can  calculate  q1  from  the  expression 


q 


f 


-K 


> 
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III-B:  Calculation  of  qls  and  gsl 

We  start  by  calculating  the  joint  probability  q,s.  The  probability  that  the  last  slot  of 
a  session  is  idle  will  be  needed  and  thus  this  quantity  is  calculated  at  first. 

Let  1^  be  a  random  variable  associated  with  the  last  slot  of  a  session  of  multiplicity 
k  that  takes  the  following  values. 

ii  _  I  0  if  the  slot  is  involved  in  a  successful  transmission 
k  - 1  1  if  the  slot  is  idle. 

From  the  operation  of  the  algorithm  it  can  be  observed  that  the  following  equations  hold. 

18  =  1  .  11=0  ,  lL  =  lU,+F2  ,  k>2.  (16) 

By  considering  the  expected  values  in  (16)  we  obtain  the  following  infinite  dimensional¬ 
ity  system  of  linear  equations 


Lb  =  1  ,  L[  =  0 

(17a) 

°°  k 

Me  =  X  X  Pf  (F2)  t>k  (<t>l)  Lk-o,+F2- 

(17b) 

F2=0  <j>,=0 


Notice  that  L*  is  the  probability  that  the  last  slot  of  a  session  of  multiplicity  k  is  idle  and 
thus  Lk  <  1  <  °°  ,  for  k<°°.  The  system  in  (17)  is  of  the  form  of  that  in  (5).  By  using  the 
same  arguments  as  those  used  in  the  calculation  of  we  can  solve  a  truncated,  up  to  J, 
version  of  (17)  and  obtain  a  very  good  approximation  of  L{..  By  averaging  the  latter  over 
all  k<J,  we  can  approximate  L\  the  probability  that  the  last  slot  of  a  session  is  idle. 

For  the  internal  IS-slots  of  the  output  process,  the  following  equations  hold. 

x{)s=0  ,  t's  =  0  (18a) 

^kS  =  ^0%F,  +  '4^0, +F2  +  l{i;|tPt  =  1  ,  k-6,  +  F2  =  1}-  (18b) 

Notice  that  the  idle  slots  which  are  the  last  of  a  session  and  are  followed  by  a  session  of 


*  v 


o 

to 
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multiplicity  1  (that  would  give  an  IS-slot),  are  not  considered  by  the  expressions  in  (18). 

By  considering  the  expected  values  in  (18),  we  obtain  an  infinite  dimensionality 
system  of  linear  equations  with  respect  to  TIS.  The  comments  that  were  made  in  the  cal¬ 
culation  of  TC1  apply  to  this  case  again  and  thus  TIS  can  be  calculated  by  solving  a  trun¬ 
cated  version  of  the  system  in  (18).  The  resulting  finite  dimensionality  system  is  of  the 
form  of  that  in  (6)  with  coefficients  akj,  0<j<J,  0<k<J,  given  by  (7)  and  constants  given 
by 

hoS  =0  ,  h\s  =  0  (19a) 

J-k  J+l-k 

hLS  =  Pf(l)  bk  (k)  £  Pf(Fi)  Lk+F]  +  Pf(0)  bk  (k-1)  £  Pf(Fi)  U_i+F,  •  (i9b) 

F,=0  F,=0 

The  average  number  of  internal  IS-slots  in  a  session,  T1S,  is  then  approximated  by 
averaging  Tks  over  all  k<J. 

By  invoking  the  strong  law  of  large  numbers  and  the  ergodic  theorem  for  stationary 
processes,  we  prove  in  Appendix  B  that  the  joint  probability  qIS  is  given  by  the  following 
expression 


’j'IS 


q“=-i-+X,e-X-  v-. 


(20) 


one 


The  joint  probability  qSI  can  be  calculated  by  following  a  procedure  similar  to  the 
developed  in  the  calculation  of  qis.  The  quantities  and  the  equations  that  are 


involved  in  this  calculation  are  given  in  Appendix  C. 


III-C:  Calculation  of  qc ,  qcc .  q"  .  qlc .  qss .  and  qsc . 

The  steady  state  probability  qc  and  the  remaining  joint  probabilities  are  calculated 


ij 


r'< 


V 

I 

>1 


from  the  following  balance  equations  of  the  probability  mass. 


In  this  section,  the  steady  state,  q  ,  and  the  joint,  q J,  probabilities  have  been  calcu 
lated,  i,  j  e  {I,  S,  C}.  The  transition  probabilities  P(i,  j)  of  the  Markov  chain  are  calcu 
lated  from  the  expression 

p(i,  j)-^r  >  i,  je{I,  S,  C}. 

q 

IV.  Performance  measure  for  the  approximation  of  the  output  process. 

In  this  section,  a  system  that  consists  of  a  central  node  which  receives  and 
retransmits  packets  originating  from  several  random  access  communication  systems,  is 
considered  (Fig.  1  ).  Each  input  stream  represents  the  output  process  from  a  multi-user 
random  access  slotted  communication  system;  that  is,  the  process  of  the  successfully 
transmitted  within  the  latter  system  packets. 

A  packet  arrival  is  declared  at  the  end  of  the  slot  in  which  the  packet  was  success¬ 
fully  transmitted.  Thus,  the  arrival  process  to  each  input  line  is  a  discrete  process.  The 
arrival  points  in  all  streams  coincide;  that  is,  the  networks  are  assumed  to  be  synchron¬ 
ized  and  all  slots  are  of  the  same  length. 


The  service  time  in  the  central  node  is  constant  and  equal  to  one,  which  is  assumed 
to  be  the  length  of  a  slot.  This  implies  that  arriving  and  departing,  from  the  central  node, 
packets  have  the  same  length.  The  first  in-first  out  (FIFO)  service  policy  is  adopted. 
More  than  one  arrivals  (from  different  input  streams)  that  occur  at  the  same  arrival  point 
are  served  in  a  randomly  chosen  order.  The  buffer  capacity  of  the  central  node  is 
assumed  to  be  infinite. 

The  objective  in  this  section  is  to  calculate  the  mean  time  that  a  packet  spends  in  the 
central  node  (waiting  time  plus  service  time),  under  the  assumption  that  the  output  pro¬ 
cess  of  a  network  is  governed  by  an  underlying  3-state  Markov  chain  and  a  correspond¬ 
ing  stationary  mapping  rule.  The  mean  time  is  used  as  a  performance  measure  of  the 
proposed  approximation  model  of  the  output  process  and  it  is  compared  with  simulation 
results  of  the  actual  system  in  the  next  section. 

The  selection  of  the  particular  system  to  measure  the  performance  of  the  proposed 
approximation,  is  due  to  the  wide  application  that  this  system  has.  For  example,  such  a 
system  appears  in  a  multi-hop  environment,  where  the  destination  of  packets  originating 
from  a  multi-user  random  access  network,  is  another  network.  Also,  such  a  system  may 
appear  in  a  single  hop  environment.  An  example  of  the  latter  can  be  a  packet  radio 
environment  in  which  more  than  one  networks  operating  in  the  same  or  neighboring 
regions  can  use  a  common  central  node  f  19]. 

A  discrete  single  server  queueing  system  with  finite  number  of  independent  input 
streams  and  per  stream  arrivals  governed  by  an  underlying  finite-state  Markov  chain,  has 
been  analyzed  in  [18].  The  system  that  is  considered  in  this  section  is  a  special  case  of 
the  general  queueing  system  in  [18].  The  state  space  of  the  Markov  chain  and  the  sta- 
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tionary  mapping  are  given  by  (1)  and  (2),  respectively. 
If  the  conditions 


N 

<  Sj(  max  ,  l^i<N  and  £  >-i  <1 

i=l 

are  satisfied,  then  the  queueing  system  is  stable.  The  average  number  of  packets  ,  Q,  in 
the  central  node  can  be  calculated  as  the  sum  of  the  solutions  of  3N  linear  equations, 
[18].  Then,  the  mean  time  that  a  packet  spends  in  the  central  node,  D,  is  given  in  con¬ 
junction  with  Little’s  formula,  by  the  following  expression. 


D  = 


(22) 


Note  that  under  stable  operation  of  the  networks  (Xj  <  Sjt  max  ,  1  <  i  <  N),  the  adopted 
mapping  rule  in  (1)  implies  that  the  input  rate  of  the  i*  stream  to  the  central  node  is 
equal  to  the  input  rate  to  the  corresponding  network. 


Let  us  denote  by  x  and  y  the  N-dimensional  vectors  that  describe  the  states  of  the  N 
Markov  chains  in  two  consecutive  time  slots,  x,  y,  e  S  =  S!x  S2x...xSN.  Let  p(j;y)  denote 
the  probability  that  there  are  j  packets  in  the  central  node  and  that  the  N-dimensional 
Markov  chain  is  in  state  y,  and  let  P(z  ;  y)  be  the  corresponding  generating  function. 
Then,  the  average  number  of  packets  in  the  system,  Q,  is  given  by  the  sum  of  the  solu¬ 
tions  of  3n  linear  equations,  [18].  These  equations  are  given  by 
N 

£  2  [2(v-l)  P  (l;x)  +  (v-l)(v-2)  P  (l;x)  +  2(v— I )  p  (0;x)|  =  0 
v=0  xeF, 

and  any  3N  -  1  from  the  following: 

N  ,  _  _  _ 

P(Ly)=Z  Z  ((v-l)P(Lx)  +  P  (l;x)  +  p(0:x)]p(x,y)  ,  y  fS  . 

v=0  *  e  Fv 


(23b) 


The  unknown  quantities  in  (23)  are  P(l;y),y  £  S;  P(l;y)  denotes  the  value  of  the 
derivative  of  P(z;y)  at  z=l.  The  set  Fv  is  given  by 

f_  -  N  ] 

Fv  =•  x  =  (xIt  ■■■xN)cS  :  Jj  ai(xi)  =  vj-. 

Since  the  input  streams  to  the  central  node  are  independent,  we  have  that 

N  N 

P(x,y)  =  riPi(*i,yi)  ,  p(0;x)  =  po  n  ^i(xi) 

i=l  i=l 

N 

P(l;x)  =  7t(x)  =  n  ^i(xi) 
i=l 

where  jtj(Xj)  and  Pi(xj,  yj)  are  the  steady  state  and  state  transition  probabilities  of  the 
Markov  chain  associated  with  the  ith  input  stream  and  po  is  the  probability  that  the  cen¬ 
tral  node  is  empty.  The  latter  is  given  by,  [18], 


Po-l-XV 

i=l 

By  solving  the  3N  linear  equations  that  are  given  by  (23)  and  summing  up  the  solu¬ 
tions,  the  average  number  of  packets  in  the  central  node,  Q,  is  obtained.  Then,  the  mean 
time  that  a  packet  spends  in  the  system  is  given  by  (22). 


V.  Results  and  conclusions 


In  this  section,  the  performance  of  the  proposed  approximation  model  of  the  output 
process  is  compared  with  the  performance  of  the  actual  system  and  that  of  some  other 
approximations.  The  mean  time  that  a  packet  spends  in  a  central  node  which  receives 
and  retransmits  packets  originating  from  N=2  and  N=3  multi-user  random  access  net¬ 
works,  is  used  as  the  performance  measure. 


w.  .  .-v. 
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We  consider  a  multi-user  random  access  slotted  communication  network  in  which 
the  binary  (C/NC)  feedback  limited  sensing  collision  resolution  algorithm,  described  in 
section  II,  is  deployed.  For  this  network,  the  output  process  is  approximated  as  described 
in  the  Introduction.  The  steady  state  and  state  transition  probabilities  of  the  Markov 
chain  are  calculated  according  to  the  procedures  developed  in  section  III. 

In  the  sequel,  we  consider  systems  of  N=2  and  N=3  networks,  as  those  described 
above,  that  feed  a  central  node.  The  mean  time  that  a  packet  spends  in  the  central  node  is 
calculated  as  it  is  described  in  the  previous  section.  The  results  (in  slots)  are  shown  in 
Tables  1  and  2,  together  with  the  results  obtained  from  the  simulation  of  the  actual  sys¬ 
tem.  It  can  be  clearly  observed  that  the  results  obtained  under  the  proposed  approxima¬ 
tion  of  the  output  process  are  very  close  to  those  obtained  from  the  simulation,  especially 
for  total  input  traffic  rate  to  the  central  node  less  than  .99. 

In  Table  3,  the  values  of  the  steady  state  and  state  transition  probabilities  that  were 
computed  from  the  procedures  developed  in  section  III,  are  compared  with  the 
corresponding  values  obtained  from  the  simulation  of  the  actual  system.  The  coin¬ 
cidence  (up  to  the  second  decimal  point  )  between  the  analytical  and  the  simulation 
results,  show  that  the  estimation  of  those  probabilities,  by  solving  truncated  systems,  is 
extremely  good.  The  difference  between  the  results  which  appears  for  small  X  and  for 
the  probabilities  that  involve  collisions,  is  due  to  the  fact  that  very  few  collisions  appear 
for  those  values  of  X  and  thus  the  simulation  results  are  not  reliable. 

If  the  output  process  of  a  network  is  approximated  by  a  process  in  which  departing 
packets  from  the  same  network  are  independent,  then  the  resulting  queueing  system  in 
the  central  node  has  been  studied  and  the  mean  time  that  a  packet  spends  in  the  central 


node,  D|,  is  given  by,  [8], 


N  N 

■ 

I  I 

n= 1  m>n 

+  N 

N 

(i-S^)L^ 

n-1 

n=l  J 

The  results  for  Dj  and  for  N=2  and  N=3  appear  in  Tables  I  and  2,  respectively. 

If  the  output  process  of  a  network  is  approximated  by  a  2-state  Markov  chain 
(arrival,  no-arrival),  then  the  resulting  queueing  system  in  the  central  node  has  been  stu¬ 
died  in  [9]  and  the  mean  time  that  a  packet  spends  in  the  central  node.  Dm  ,  is  given  by 


Dm  - 


Z  Z  KK  i  +  + 

n=l  m>n  1  in  1  im 


N  N 

n=l  n=l 


where  ya  =  P(l/1)  -  P(l/0)  and  where  P(l/1)  and  P(l/0)  denote  conditional  probability  of 
arrival  given  arrival  in  the  previous  slot,  and  arrival  given  no-arrival  in  the  previous  slot, 
respectively.  Clearly,  DM  can  be  also  obtained  by  solving  2N  equations,  according  to  the 
procedures  described  in  the  previous  section.  The  results  for  DM  and  for  N=2  and  N=3 
appear  in  Tables  1  and  2,  respectively. 

By  examining  the  results  in  Tables  1  and  2,  we  conclude  that  when  the  total  input 
traffic  to  the  central  node  is  less  than  .99  and  the  input  traffic  to  each  network  is  substan¬ 
tial  (  >  .25),  then  the  proposed  approximation  is  the  best  among  all  considered  in  this  sec¬ 
tion.  The  proposed  approximation  performs  also  well  for  light  traffic.  The  2-state  Mar¬ 
kov  approximation  seems  to  perform  better  for  light  input  traffic  to  each  network.  In  this 
range  of  input  traffic,  the  simulation  results  are  not  very  reliable  since  arrivals  and  espe- 
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daily  collisions  are  very  rare.  Furthermore,  the  differences  among  the  results  is  so  small 
that  conclusions  is  not  easy  to  be  drawn.  On  the  other  hand,  for  substantial  input  traffic 
to  each  network,  the  2-state  Markov  model  fails  completely.  This  was  expected  since  the 
latter  model  fails  to  identify  idle  slots  from  collided  and  there  is  a  substantial  number  of 
the  latter  in  the  output  process. 
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Table  2.(N=3) 


.10 _ 30 _ 30 _ .33 

Anal-Sim  Anal-Sim  Anal-Sim  Anal-Sim 

.100-. 100  .200-300  [  300-.  300  '.330-330 

.887-. 887  .730-. 731  .460-.461  327-329 

.013-.013  .070-  069  3.40-.239  “  343-341 

.215-. 212  .184-. 184  .156-. 155  .148-.  147 

,463-.466  .  424-. 424  384-385  372-373 

322-323  .392-392  .460-.460  .480-.480 

.845-. 846  .665-. 668  .408-.410  ,294-.297 

.138-. 137  .256-. 252  359-358  390-388 

.017-.017  t~.079-.079  .23.3-.232  316-315 

.902-. 902  .800-. 799  .653.653  .548-.550 

.090-. 090  .164-.  164  .217-  219  I  .225-325 

.007-. 007  j  .037-.036  ,130-.128~  .227-325 

*  .9899-. 9899  **  .0001 -.0001 

***  .00016-00007  ****  .000075-. 00007 

Table  3. 


From  the  description  of  the  algorithm,  the  following  equations  can  be  written. 


To  =  1  ,  T1!  =  0 

Tk  =Toi+F,  +Tk-0,+F2  ,  k>2 
and 

XqS  =  0  ,  tP=0 

TkS  =  1{6-+F,  =  1}  +  T^+F,  +  xk-0i+F2  > 

where  the  random  variables  involved  have  been  defined  before.  By  considering  the 
expectations  in  the  above  equations  and  by  truncating  the  resulting  system  of  linear  equa¬ 
tions,  we  obtain  a  finite  systems  of  the  form  of  (6)  with  coefficients  gi  /en  by  (7)  and  con¬ 
stants 

ho  =  1  ,  h[=0  ,  k>l 
and 

hP=0  ,  hfs=0  ,  hP=Pf(0)bk  (1)  +  Pf(l)bk(0)  ,  k>2. 

Appendix  B 

In  this  Appendix  we  prove  equation  (20).  Let  us  define  the  following  random  vari¬ 
ables. 

j is  I  1  if  the  n^slot  of  a  session  is  an  IS— slot 
n  1 0  otherwise 

NV 

j|,,s  =  £  I„s  :  Number  of  IS— slots  of  the  vIh  session 

n=l 

Nv  :  length  of  the  Vth  session  (in  slots) 


e  is _ J  1  if  the  vlh  session  has  a  non-imemal  IS— slot 

~  1  0  otherwise 


Notice  that  is  associated  with  the  last  slot  of  the  Vth  session.  A  non-internal  slot  is 
assigned  to  the  session  which  its  first  slot  belongs  to. 

The  joint  probability  qIS  can  be  calculated,  by  using  the  law  of  large  numbers,  from 
the  following  expression. 

M 


I  JIS 


IS  i'  Vs! 
q  =  lim 


M 


M 


INv 

V=1 


The  above  expression  can  be  written  as 
M 


IS, 


qls  =  lim 

M— >~> 


v=l 


M 

Z  Nv 

V=1 


+  lim 


M 

141s 

v=l 


M 


M 


I  Nv 
v=l 


M 


£2* 


M 


=  lim 


V=1 


M 


14 


V=1 


M_>"  15-  s  N. 

M  V=1 


+  lim  .. 
M— >°°  1  M 


M  5N” 


Clearly,  the  random  variables  Nv,  tyS  ,  v>l,  are  independent  and  identically  distributed 
with  mean  value  L  <  °°  and  TIS  <  «»  (since  TIS  <  L  ) ,  respectively.  Thus,  the  strong  law 
of  large  numbers  asserts  that 

M 


q,S  = 


•pIS  M 


1  M  1C 

lim  77  I  4v 

M 


+ 


v=  1 


(Bl) 


L 


The  random  variables  ,  v>  1,  are  not  independent  but  {^[?}v>i  is  a  stationary  pro¬ 
cess  and  £yS  has  expected  value  given  by 


26 


E{^lS }  =  P{ (I— slot  ,  k  =1)  /  last  slot  of  a  session} 

where  k  denotes  the  multiplicity  of  the  next  session.  Since  the  two  events  are  indepen¬ 
dent  we  can  write  that 

E{^s}  =  Ll  <=»  .  (B2) 

The  second  term  of  the  above  product  is  the  probability  of  having  a  session  of  multipli¬ 
city  1.  By  applying  the  ergodic  theorem  for  stationary  processess  to  (B 1 ),  [17],  and  by 
considering  (B2),  we  obtain  the  expression  in  (20). 


?endix  C 


In  this  Appendix  we  present  the  procedure  to  calculate  the  joint  probability  qSI.  Let 
lk,  Lk  and  Ls  be  defined  as  the  quantities  lk,  L^  andL1,  respectively,  by  replacing  the 
term  idle  by  the  term  success.  Since  the  last  slot  of  a  session  is  either  idle  or  involved  in 
a  successful  transmission  we  have  that 


L|  -  1  - 
Ls  =  1  -V. 

The  equations  that  correspond  to  those  in  (18)  are  given  by 

, t§'=o 

Tk'  =  toJ+F,  +  tk-0,+F2  +  l{i;,.Fl  --  1  ,  k-O.+Fj  =  0}  -  k  >  2 
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